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Abstract. By means of population synthesis we find that 
, the expected Galactic number of black holes with massive 
helium star companions is ~ 100 and depends on an as- 
sumed threshold for M pro _BH- The overwhelming major- 
ity of these systems has orbital periods in excess of 10 hr, 
with a maximum at ~ 100 hr, while under the Illarionov & 
Sunyaev (1975) disk formation criteria for accretion from 
, the strong stellar wind of Wolf-Rayet star disk accretion is 
possible only for orbital periods below ~ 10 hr. However, 
, the number of such short-period systems is vanishingly 
' small. If the accretor in Cyg X-3 is a 10 M© black hole, 
then the accretion rate will be super-Eddington. Super- 
Eddington accretion may be responsible for the formation 
, of jets in Cyg X-3 and may also support an X-ray lumi- 
nosity as high as ~ 10 39 erg s — 1 . From the orbital pe- 
riod distribution for neutron stars with massive helium 
' companions we find that if during the common envelope 
phase a neutron star accretes at AfEdd and spins-up to 
' the equilibrium period, then in most systems it acts as 
, a "propeller" and accretion from the WR star wind is im- 
possible. For the model with two massive helium stars as 
, an immediate progenitor of Cyg X-3, requirement of acco- 
• modation of two WR stars in the post-common-envelope 
, orbit combined with severe mass loss by them prevents 
formation of BH+WR systems with orbital periods less 
than several days. 



Key words: binaries: close-binaries: accretion-stars: 
Wolf-Rayet 

1. Introduction 

At present, three black hole candidates with massive early- 
type companions are known (Cyg X-l, LMC X-3 and LMC 
X-l). About half a dozen of black hole candidates with 
low-mass companion stars have been identified as the so- 
called X-ray transients (Tanaka &; Lewin 1995; Tanaka & 
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Shibazaki 1996). Yet another candidate may be the puz- 
zling 4.8 hr orbital period X-ray binary Cyg X-3 (Giacconi 
et al. 1967). Its spectrum is hard, with a tail extending to 
ultra-high energy 7-rays (Lloyd-Evans et al. 1983). Cyg 
X-3 shows an asymmetric sinusoidal modulation in the 
low-energy X-ray emission (Mason et al. 1976). The same 
modulation is present in the infrared light curve (Beck- 
lin et al. 1973; Mason et al. 1986). Infrared spectroscopy 
has recently shown that the companion is a Wolf-Rayet 
star of the WN7 subtype (van Kerkwijk et al. 1992, 1996; 
van Kerkwijk 1993). Cyg X-3 has huge radio outbursts 
with evidence for jetlike emission expanding at ~ 1/3 of 
the light velocity (Gregory et al. 1972; Geldzahler et al. 
1983). 

The orbital period of Cyg X-3 is increasing on a 
time scale of 850000 years (Kitamoto et al. 1995). 
With this P/P it is possible to estimate the stellar 
wind mass loss rate as Md yn = M tot (P /2P) ps 6 x 
10- 6 (M to t/10M Q ) Mq yr" 1 , where M tot = M c + M W r, 
M c and Mwr are the masses of the compact object and 
the WR star, respectively. 

Another estimate of the mass loss rate can be made 
using the infrared flux Fg. Van Kerkwijk et al. (1996) 
estimated Mg using a simple relation Mg /v w oc Fl'\ For 
the wind velocity of 1450± 150 km s _1 and Fg determined 
from the May 1992 infrared light curve they find Mg ~ 
1.2 x 10~ 4 M yr _1 , an order of magnitude larger than 
the dynamical estimate Mdyn (for Mwr =10 M©, M c = 
1.4M©). But as was already discussed by van Kerkwijk 
et al. (1996), Mg ~ 10~ 4 M yr" 1 that is inferred from 
the infrared flux should be regarded as highly uncertain. 
Due to dumpiness of the wind, the true mass loss rate 
may be lower than Mg by a factor 2-3 and for Cyg X-3 
the difference may be even larger (see also Cherepashchuk 
& Moffat (1994) and references therein for a discussion of 
the effect of clumps). 

Cyg X-3 is a strong X-ray source with mean X- 
ray flux F x (2 - 20KeV) » 6.8 x 10~ 9 erg s^cm- 2 , i.e. 
L x » 1.2 x 10 38 (d/12 Kpc) 2 erg s^ 1 (Kitamoto et al. 
1995). Cherepashchuk & Moffat(1994) argue that since 
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the X-rays is observed in the IR range, the true intrinsic 
X-ray luminosity of the accreting relativistic object in Cyg 
X-3 must be considerably higher than the observed mean 
value of the hard X-ray luminosity Lx ~ 10 38 erg s _1 . 
In their opinion, "this fact favours the presence of an ac- 
creting black hole, as opposed to a neutron star, in Cyg 
X-3". 

Recently, Schmutz et al. (1996) observed time- 
variations in the profiles of several infrared emission lines 
from Cyg X-3. They concluded that the variations are due 
to the orbital motion of the WR star and derive a mass 
function for Cyg X-3 of 2.3 M@. Assuming reasonable val- 
ues for the mass of WR star and the inclination angle, 
they obtained a range of masses 7-40 Mq for the compact 
object, with the most likely value of 17 Mq, from which 
they concluded that the compact component of Cyg X-3 
is a black hole. 

Van den Heuvel & De Loore (1973) were the first to 
suggest that Cyg X-3 may be a system in a post massive X- 
ray binary (MXRB) stage and may harbour a compact ob- 
ject (neutron star or black hole) and a helium star of sev- 
eral solar masses. Tutukov & Yungelson (1973a) indepen- 
dently suggested that the optical components of MXRB 
will in the later evolutionary stage become WR stars pro- 
vided that they are massive enough. 

In the present paper we attempt to model the Galac- 
tic population of black- hole + massive helium star systems 
and to show that the very short orbital period makes Cyg 
X-3 an extremely rare example of such systems detectable 
by X-ray observations. Also, we discuss the circumstances 
which prevent discovery of a large Galactic population 
of neutron star+helium star systems, which are probably 
quite numerous. 

2. Cyg X-3 - a unique system with a black hole 
accretor? 

2.1. Formation scenarios for BH + WR binaries 

The calculations for the present study were made by 
means of the population synthesis code which was pre- 
viously applied by one of the authors (LRY) for the mod- 
elling of different components of the binary star popula- 
tion of the Galaxy, among them binary Wolf-Rayet stars 
(Yungelson & Tutukov 1991; Vrancken et al. 1991), neu- 
tron stars (Tutukov & Yungelson 1993a,b) and high-mass 
X-ray sources (Iben et al. 1995). The basic assumptions 
important for the present study are the following. 

The initial mass function of the primaries in binary 
stars is a power-law dN oc M^ 2 ' b dM\. The initial distri- 
bution of binaries over separation a is flat in log a for < 
log(a/i? Q ) < 6. The orbits of close binaries are initially 
circular. The initial distribution of binaries over mass ra- 
tios of components go = M<z/M\ is flat for < Qq < 1. 
Discussion of the effects of the assumed distributions over 



a and qo upon the population of high-mass binaries in the 
Galaxy may be found elsewhere (e.g. Portegies Zwart & 
Verbunt 1995; Portegies Zwart & Yungelson 1997). They 
appear to be of no serious significance to the (mostly) 
qualitative conclusions of the present paper. 

Normalization of the Galactic stellar birthrate function 
corresponds to the formation of one binary system with 
the initial mass of the primary greater than 0.8 Mq per 
year. The rate of binarity of stars is assumed to be equal to 
100 %.□ With this normalization and the assumed rate of 
binarity our model gives a combined rate of SN Ib/c and 
SN II ~ 0.02 yr" 1 , consistent with the sum of expected 
rates of Galactic SN Ib/c and SN II (Cappellaro et al. 
1997). 

Below, we use the following notation: MS - main- 
sequence star, RLOF - stage of the Roche lobe over- 
flow (stable or with formation of a common envelope), He 

- helium star remnant of a star which experienced a mass 
loss (in fact, a WR star) , SN - supernova explosion, BH 

- black hole. Then all basic scenarios for the formation of 
black-hole + helium star binaries may be summarized as: 

A. MSi + MS 2 -> RLOFi -> Hei + MS 2 -> SNi + MS 2 
-> BHi + MS 2 -> RLOF 2 -> BHi + He 2 



and 

B. MSi+MS 2 - 
-> BHi + MS 2 



Hei + MS 2 - 
> BHi + Hc 2 . 



SNi + MS 2 



Scenario A occurs in the systems in which stellar wind 
mass loss does not prevent significant expansion of stars 
and RLOF. In the opposite case scenario B occurs. In 
fact, in scenario B components of binaries evolve inde- 
pendently like single stars. In absence of hydrodynamic 
calculations of the mass exchange, the upper mass limit 
of stars which may experience classical RLOF is, in fact, 
one of the free parameters of the problem (for discussion 
of this issue see, e.g. Sybesma 1986; Moffat 1995). Con- 
ventional computations of stellar evolution with stellar 
wind mass loss rates suggested by observations (e.g. Mas- 
sevich et al. 1979; Sybesma 1986; Vanbeveren 1991) show 
that scenario A occurs in systems with the initial mass 
of the primary component below 35 - 50 Mq. However, 
depending on the mass accumulated in the RLOFi event, 
the secondary in scenario A may avoid RLOF and, vice 
versa, the secondary in scenario B may encounter RLOF 
if its mass is sufficiently low. In the present study we as- 
sume that the border between scenarios A and B is at 
Mi =50 Mq. Systems which evolve under the influence 
of stellar wind mass loss have larger orbital separations in 
the BHi + He 2 stage than systems which evolved through 
nonconservative RLOF. Decrease of the abovementioned 



For another birthrate normalization and assumptions on 
the rate of binarity all numbers of binaries given below have 
to be rescaled. 
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border below Mi = 50 Mq would make the conditions for 
the formation of Cyg X-3-like systems even less favourable 
than is found below (see Fig. 1). 

The mass of the helium star immediately after the 
end of mass exchange episode is given by the relation 
Mh c = 0.066Mj^g 4 , which represents a fit to the evolu- 
tionary calculations by Tutukov & Yungelson (1973b) and 
Iben & Tutukov (1985). 

To estimate the variation in the separation of compo- 
nents in the common envelope stage we applied the equa- 
tion suggested by Tutukov & Yungelson (1979): 

where subscripts and / refer to the initial and final 
states, respectively, and a ce is the parameter which mea- 
sures the ratio of the binding energy of the common en- 
velope and the change of the orbital energy of the binary 
(see e.g. Iben & Livio 1993; Livio 1996 for discussion). 
This equation implies that the spiral-in starts when the 
envelope surrounding the two cores expands to ~ 2ao and 
mass is actually ejected from the gravitational potential of 
both components. Computations were performed for the 
common envelope parameter a ce — 1, which gave reason- 
able results in our previous attempts to model the popu- 
lation of binary stars. When compared to the equation for 
(if /do derived by Webbink (1984), which is often applied 
for population synthesis studies, a ce = 1 in Eq. (1) corre- 
sponds to a ce k 4 in Webbink's equation. For the latter 
it formally means that most of the energy for expulsion 
of the common envelope comes from sources other than 
the orbital energy of the binary. However, there are clear 
indications that with a ce ^3 — 4 (in Webbink's formula- 
tion) one would obtain more reasonable agreement with 
observations for several classes of descendants of massive 
binaries than with lower values of a ce (Portegies Zwart & 
Yungelson 1997). 

2.2. Masses of the Cyg X-3 components 

The actual mass range of progenitors of black holes both 
in single stars and in binaries is, in fact, unknown. The 
main factors contributing to this are the uncertain mass 
loss rates of hydrogen- and helium-rich stars, poor un- 
derstanding of mixing processes in stellar interiors and 
of processes during a supernova explosion and a poorly 
known equation of state (see e.g. Woosley ct al. 1995; 
Brown et al. 1996 and references therein). The "standard" 
assumption is that binary components more massive than 
40 Mq form black holes (van den Heuvel & Habets 1984). 
Recent calculations of Brown et al. seem to confirm this 
limit, while according to Woosley et al. (1995), who as- 
sume an extremely high mass-loss rate for helium stars, 
this limit may be as high as ~ 60 Mq . Another impli- 
cation of the Woosley et al. results is convergence of pre- 
supernovae masses to 3-4 Mq . Actually, significant spread 



of masses in observed candidate black holes in soft X-ray 
transients (4-16 Mq) suggests that there may be factors 
other than the initial mass which determine the fate of 
a star (Ergma & van den Heuvel 1997). Under the given 
circumstances, we treated the masses of black holes and 
their progenitors as an additional model parameter. Com- 
putations were performed on the assumption that black 
holes result from objects which have a mass mass exceed- 
ing 7 or 10 Mq at the end of the helium star stage. Under 
the initial-final mass relation given in Sect. 2.1 and us- 
ing Langer's (1989) mass loss rates for helium stars, these 
limits correspond to Mzams ~ 30 or 50 Mq . We assumed 
that no mass is ejected during a massive helium star rem- 
nant collapse into BH (Burrows 1987; Woosley & Weaver 
1995). 

The mass of the WR component of Cyg X-3 is like- 
wise uncertain. Cherepashchuk & Moffat (1994) suggest 
the range 10-50 Mq based on the spread of estimates 
of masses of actually observed WN7 stars. One can use 
the estimated mass loss rate of the stellar wind to get 
an idea about the WR star mass in Cyg X-3. From 
the observations for five binary Wolf-Rayet stars, Abbott 
et al. (1986) derive M w oc M^\. Langer (1989) sug- 
gests, from the fits of stellar models to the observational 
data, that in the WR star core helium burning phase 
M w w (1 - 6) x 10- 8 (M W r/M q ) 2 - 5 M q yr" 1 . If M w ~ 
l/10M ff w 1.2 xlO" 5 M yr" 1 (van Kerkwijk et al. 1996) 
and M w w Md yn , then M to t ~ 20M Q . A direct application 
of Langer's formulae gives 8 ^ Mwr/M q 17. 

The rate of accretion onto a black hole may be esti- 
mated as 



47ra 2 



where r a 



2GM] 



BH 



(2) 



Then, according to Kepler's third law 



Mace « 0.14. 



A/f 2 4 



,4 

1000 



-M WR )3M dyn , (3) 



where all masses are in Mq, Ph r is the orbital period in 
hours, and wiooo is the wind velocity in 1000 km s _1 . For 
M BH = 10M Q , t> 1000 = 1.5, P hr = 4.8, M WR = 10M Q 
one gets M acc w 0.05 x M dyn w 5.7 x 10~ 7 M yr" 1 . 
The Eddington limit for the accreting black hole (helium 
accretion and R = 3r g ) is 2.3 x 1O~ 8 (Mbh/M ) M yr" 1 
and for M B h = 10M Q it is M E dd « 2.3 x 10~ 7 M Q yr _1 . 
Thus, in the latter case the accretion rate exceeds the 
Eddington limit. It is easy to find that if the accreting 
object is a neutron star, then 

M acc (NS) = f^Y ( M ™ + M ™\ % 

accV ' \MbhJ \M wr + MbhJ 1 ' K ' 

For the same mass of the Wolf-Rayet star (10 Mq), but for 
a neutron star or low-mass black hole accretor (1.5 Mq, 
Brown 1995), it is simple to estimate by Eq. (4) that M acc 
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is below the Eddington rate. However, if the accretion is 
super-Eddington, it may be responsible for the formation 
of jets in Cyg X-3. A high accretion rate may support an 
X-ray luminosity as high as ~ 10 39 erg s _1 which is con- 
sistent with Cherepashchuk & Moffat's (1994) estimate. 

No kicks were imparted to the nascent BH. Even in 
the case of neutron stars the problem of natal kicks is still 
unsolved: observational data may be explained both on the 
assumption of the absence of kicks (Ibcn & Tutukov 1996) 
and of their presence (van den Heuvel and van Paradijs 
1997). As an argument for absence of kicks for BH one 
may consider the fact that black hole candidates in low- 
mass X-ray binaries (LMXB) share the radial velocities of 
their local rest frames in the Galaxy (with the exception 
of Nova Sco 1994). The dispersion of their distances to 
the Galactic plane is smaller than for LMXB with neutron 
stars by more than a factor 2. This indicates that any kick 
velocities that the black hole systems may have received 
at their formation are considerably smaller than those of 
neutron star systems (White & van Paradijs 1996). 

2.3. Model population of BH + WR and NS + WR bina- 
ries 

Figure 1 shows the distribution of BH + WR binaries over 
orbital periods and masses of WR stars for two assump- 
tions on the progenitor mass of the BH. The minimum 
mass of the WR star is taken equal to 5M . The difference 
between the two distributions may be due to the fact that 
for M pre _BH > 30 M Q both scenarios A and B operate, 
while for M pre _BH > 50 M Q only scenario B is effective, 
i.e. in the first case some of pre-WR stars were able to 
accumulate some mass via stable RLOF. The predicted 
numbers of BH + WR stars in the two cases are ~ 300 
and ~ 100, respectively. These numbers can be considered 
only as giving an order of magnitude estimate on the po- 
tential number of the Galactic BH+WR systems. If natal 
kicks or mass ejection accompany the formation of BH, 
they may reduce the number of these systems, but they 
cannot be overly significant for the most massive binaries 
we are dealing with. 

More important, in our opinion, is the fact that the 
overwhelming majority of these systems have orbital peri- 
ods longer than 10 hours, with a maximum of ~ 100 hours. 
This may be helpful in to understanding why Cyg X-3 is 
unique. On the basis of the seminal Illarionov & Sunyaev's 
(1975) paper "Why the Number of Galactic X-Ray Stars 
Is So Small?" one can make the following estimate. Illari- 
onov & Sunyaev have shown that in the case of accretion 
of stellar wind matter in a detached binary system the 
specific angular momentum of the matter captured by the 
relativistic star is typically small. Therefore, no accretion 
disk is formed around the relativistic star. Consequently, 
very special conditions are required for a black hole in a 
detached binary system to be a strong X-ray source. A 
disk may form if the specific angular momentum of ac- 
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Fig. 1. Normalized number-density distribution of BH+WR 
systems over orbital periods and masses of WR components 
for two assumptions on the minimum mass of the progeni- 
tors of black holes - 30 Mq (upper panel) and 50 Mq (lower 
panel). Dashed line labelled by 1.0 gives the upper limit of 
orbital periods for which disk accretion onto BH is possi- 
ble if Mbh = WMq and viooo = 1-0. Maximum of the 
gray-scale is the same for both panels and corresponds to 
aiogtiaiogP ~ 274, where N is the number of systems. 



creting matter exceeds the specific angular momentum of 
the particle in the last inner stable orbit with R — 3r g 
which is equal to Q m i n = y/3r e c. Thus, the disk formation 
criterion is 



-nRi > V3. 



■rvc, 



(5) 



where R\ is the capture radius of stellar wind matter by 
the relativistic star. This condition is fulfilled when the 



E. Ergma & L. Yungelson: Cyg X-3 



5 



NS+WR Binaries 



NS+WR Binaries 



Scale 



M pre _ BH > 30M o 



4 - 




1.0 
1.5 



0.6 



1.2 



1 .4 



1 .0 



Mpre-BH * 50M Q 



2 - 



0.6 




0.8 1 1.2 

log (M WR /M Q ) 



. 1.0 
- 1.5 



1 .4 



= 1.0 



M pre _ BH > 3OM 



1 500km/s 



4 - 




_ -12 
„ -1 1 
- - '10 



0.6 



1.2 



1 .0 



Mprs-BH ^ 50M G 

v w = 1500km/s 



2 - 




0.6 



0.8 1 1.2 

log (M WR /M ) 



1.4 



-12 
-1 1 

10 



1.4 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




Fig. 2. The same as in Fig. 1, but for NS+WR systems. In panels (a) and (b) dashed lines labelled by 1.0 and 1.5 give the upper 
limits of orbital periods for which "propeller" effect does not prevent disk accretion onto NS if t>iooo = 1-0 and 1.5 respectively. In 
panels (c) and (d) dashed lines labelled 10, 11, and 12 give the upper limits of the the orbital periods for which spin deceleration 
timescale of a neutron star is shorter than the helium-burning time of its companion if magnetic field strength is B = 10 10 , 10 11 , 
and 10 12 Gs, respectively. Maximum of the gray-scale is the same for all panels and corresponds to 81og ^aiog p = H60. 



period of binary system P or b < 4:.8(Mbh/ Mq)viq Q0 5 2 nr ; 
where 6 <~ 1 is a dimensionless parameter. Typical Wolf- 
Rayet stellar wind velocities range from 1000 km s _1 to 
4000 km s _1 (Conti 1988). As a result, disk accretion is 
possible only for systems having a very short orbital pe- 
riod (Fig. 1). Curiously enough, if Mbh ~ 5M and wind 
velocity <~ 1500 km s _1 for Cyg X-3, the disk accretion 
condition is satisfied when P or b ~ 4.8 hr! (A similar con- 
clusion was reached by Ibcn et al. 1995). 

A Wolf-Rayet binary with P or & = 432 and sus- 
pected BH-companion HD 197406 (Cherepashchuk 1991; 
Marchenko et al. 1996) may be a representative of the 



population of BH+WR binaries which do not show up as 
X-ray sources. 

The population of BH + WR stars can be com- 
pared with the population of neutron stars accompanied 
by WR stars (henceforth, NS + WR). Figure 2 shows 
logMwR — logP plots for NS + WR stars similar to Fig. 
1. The respective numbers of systems are ~ 570 and <~ 700 
for Mp r e_NS < 30 and 50 Mq. These numbers provide up- 
per limits for this possible population because no natal 
kicks were imparted to neutron stars. It is also worthwhile 
to note that natal kicks tend to influence more strongly the 
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eccentricities of orbits of post-supernovae systems than 
their orbital periods. 

One may ask why NS + WR systems are not observed. 
We would like to call attention to the following. If the WR 
component is formed via RLOF one can expect that a 
common envelope forms because of the high mass ratio of 
components at the instant of the RLOF by the progenitor 
of the WR star. As a rule, the common envelope stage 
lasts for <~ 10 3 — 10 4 yrs. During this stage the NS may 
accrete at the Eddington's rate of 1.5 x 10~ 8 M Q yr -1 . 
Accretion spins up the neutron star (Pringlc & Rccs 1972; 
Davidson & Ostriker 1973) to the equilibrium period 

P ro t=P eq ~4Bf 1 M^M-i s, (6) 

where B„ = B/10 11 G, M n = M/1.5 x 1CT 11 M Q yr" 1 , 
Mns is the neutron star mass in Mq . The amount of mat- 
ter required for the spin-up to P eq is 

A(M/M )«O.l(Peq/l-5ms)-*. (7) 

For B\\ — 1, 10 and M\\ = 1000, the equilibrium peri- 
ods are 0.16 and 1.17 s, respectively. The amount of matter 
required for accretion is for these two cases AM = 1.9 x 
10~ 4 Mq and 1.4 x 10~ 5 M . These amounts may be eas- 
ily accreted during the common envelope stage. After the 
common envelope is dispersed, the NS appears to be im- 
mersed in the strong wind of WR star. Its initial rotational 
period is determined by Eq. (6) for accretion at MEdd- 
Following Illarionov & Sunyaev (1976) (see also Lipunov 
1982), accretion onto the surface of a rotating neutron star 
is possible only if its rotational period is P rot £ P{, q , where 
P c ' q is now an equilibrium period determined by the rate of 
accretion from the WR star wind. Then, from Eqs. (2) and 
(6), for M W r = 4 x 10- 8 M^ R M yr" 1 (Langer 1989), 
Msdd = 1.5 x 10~ 8 M Q yr -1 , it follows, that accretion is 
possible only for systems with 

P olh < 7.3B^M^ s M$ R M t -jP^ l hr (8) 
or 

P orb < 0.46Mj s Mj M# R «r 3 00 hr, (9) 

where M to t = Mns + Mwr- For a higher P or b a neutron 
star would act as a "propeller". Note, that the critical 
value of P rb docs not depend on B [in eq.(9)] after elimi- 
nation of the equilibrium period of rotation in the common 
envelope. For t>iooo = 1 an d 1.5 the limiting P or b is plotted 
in Fig. 2 [panels (a) and (b)]. It is evident that for the over- 
whelming majority of NS + WR systems the "propeller" 
effect may prevent accretion, hence, X-rays emission. This 
result does not depend on the value of a ce assumed in the 
population synthesis calculations, because the lower limit 
of the orbital periods of systems under consideration is 
determined by the requirement of accomodation of the 
WR star in the post-common-envelope orbit. Application 



of other prescriptions for mass and momentum loss in the 
evolution of close binaries cannot change this result ci- 
ther, because the only difference which can be expected is 
a different distribution over orbital periods with the same 
minimum period [see Portegies Zwart & Verbunt (1995) 
for a discussion of the influence of different population 
synthesis "recipes" ] . 

A neutron star in the "propeller" stage would experi- 
ence spin deceleration. The characteristic time of a neu- 
tron star spin-down is (Illarionov & Sunyaev 1975) 

i A «4x lO^rM^M-^fooo yr (10) 
or 

_3 12 55 11 22 45 

t A « 3.8 x 10 6 P n ? M NS 7 M w -M to «P o -< V> Yr- (H) 

This deceleration time can be compared with a helium 
burning time in the core of a WR star for which the fit 
to the Paczyhski (1971) and Ibcn & Tutukov (1985) data 
gives: 

logt Hc ~ 7.15-3.7?/ + 2. 237/ 1 - 37 yr, (12) 

where y — log(MwR/M ). The spin of the neutron star 
decelerates enough to allow accretion if orbital period of 
the system is 

21 9 15 3 1 35 

Porb Z OXigBg M^ R M t >-^ hr, (13) 

where tn e is in units of 10 6 yr. This critical period is plot- 
ted in Fig. 2 [panels (c) and (d)] for WR wind velocity 1500 
km s _1 and three values of the strength of the magnetic 
field which fit the range expected for neutron stars of the 
age of several Myr, typical of NS + WR systems. Figure 2 
shows that even if deceleration is taken into account, the 
propeller effect excludes accretion for the overwhelming 
majority of NS +WR binaries. 

If a WR star forms due to wind mass loss, one can 
expect a much lower accretion rate onto NS than in the 
common envelope, hence, a longer equilibrium period. For 
the same wind velocities, conditions for disk formation 
may become more favourable. However, as Fig. 2 clearly 
shows, descendants of very massive stars probably com- 
prise a small fraction of all WR stars in NS + WR binaries: 
e.g. \ogMwnl Mq = 1.2 under our formalism corresponds 
to a main-sequence mass of ~ 35M . 

2.4- Future evolution of the Cyg X-3 system 

Accepting that Cyg X-3 contains a 10 M WR star and 
a 10 Mq black hole it is possible to speculate about pos- 
sible endpoints of its evolution. A substantial amount of 
mass may be lost in the WR evolutionary stage. This mat- 
ter leaves the system in an isotropic, high-velocity wind, 
carrying away specific angular momentum of the primary. 
Then 

£ = Mm (i4) 

a Mbh + Mwr 
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For a 10 M Q helium star, the evolutionary timescale is 
<~ 5.5 x 10 5 yr (Langer 1989). At the end of the evolu- 
tion of the present WR star, system will have ft « 8 
hr, MwR,f ~ 5M Q . If the outcome of evolution is a col- 
lapse into a black hole, then the system's further evolution 
will proceed due to gravitational wave radiation and in 
~4x 10 8 yrs time the system will merge. If a Wolf-Rayet 
star forms a neutron star and expels its envelope, then, 
depending on the magnitude and direction of the natal 
kick, its semi-major axis may become so large that the 
neutron star will not merge with its black hole companion 
in a Hubble time. 



3. Discussion and conclusion 

Our study, as one could expect from the previous stud- 
ies of a similar kind, suggests that there may exist a sig- 
nificant (<~ 100 objects) Galactic population of massive 
(<~ 10M Q ) black holes and neutron stars accompanied by 
Wolf-Rayet stars. Yet, the only example of such systems 
observed to-date is Cyg X-3 which is suspected to har- 
bour a black hole. Rather rough estimates made in Sec. 
2 suggest a combination of reasons for the low propor- 
tion of such observed systems. First, both BH+WR and 
NS+WR systems predominantly form with their orbital 
periods exceeding several days. Second, for disk accretion 
from WR star wind to occur, the angular momentum of 
the wind matter has to be high. Hence, the formation of 
accretion disks is possible only in BH+WR systems with 
relatively short orbital periods. However, the number of 
such systems is vanishingly small (Fig. 1). Cyg X-3 may 
be a unique example of such a system. 

The formation of accretion disks in NS+WR systems 
is excluded by their specific evolutionary history. Namely, 
the spin-up of neutron stars due to Eddington rate accre- 
tion in common envelopes which accompany the birth of 
the overwhelming majority of WR stars suffices to prevent 
accretion via the action of the "propeller" mechanism, if 
in the subsequent stage, the wind velocity of the WR star 
exceeds ~ 1000 km s _1 (Fig. 2), the actual lower limit for 
WR star winds. The lower limit of the model orbital pe- 
riods of BH/NS+WR systems depends solely on the radii 
accepted for helium stars and does not depend on other 
parameters of population synthesis models, such as the 
common envelope parameter. 

Recently, the possibility that Cyg X-3 may harbour 
a black hole as the compact object has been discussed 
by Brown(1995). Having considered the masses of compo- 
nents in this system, Brown (1995) suggests that the pro- 
genitor of the present compact object had Mzams ~ 35 
Mq. The immediate progenitor of Cyg X-3 in Brown's 
model contains two massive WR stars. However, the re- 
quirement of accomodation of two WR stars in the post- 
common-envelopc orbit combined with severe mass loss 
by helium stars prevents, in this model, formation of 



BH+WR systems with orbital periods less than several 
days. 

It is also suggested that low-mass black holes in close 
binaries may originate as a result of hypercritical accretion 
onto neutron stars in common envelopes (Zel'dovich et al. 
1972; Chevalier 1993; Brown 1995; Fryer et al. 1996). How- 
ever, it is unclear whether these black holes may survive 
in the common envelopes and avoid swallowing of whole 
envelope of the companion. 

All low-mass X-ray binaries with suspected black hole 
components (with a possible exception of GX 339-4, Mak- 
ishima et al. 1986) are soft X-ray transients (SXT), while 
SXT are absent among massive X-ray binaries. Since Cyg 
X-3 has a short orbital period, similar to the periods 
of SXT's, one may wonder why Cyg X-3 is a persis- 
tent source. Recently, van Paradijs (1996) and King et 
al. (1996, 1997) have discussed the conditions necessary 
for disk instability in low-mass X-ray binaries in the pres- 
ence of X-ray irradiation. King et al. (1997) have realized 
that irradiation is much weaker if the accreting object is a 
black hole rather than a neutron star, since black holes do 
not have hard surfaces and cannot act as point-like central 
sources. According to King et al. , a system is persistent 
in X-rays if 

M > M^ t « 2.86 x IO-^MbJm^pI M q yr" 1 , (15) 

where M2 is the mass of the companion to the BH. For 

Mbh = M 2 = 1OM , P hr = 4.8, MZt « 10~ 9 M yr" 1 , 
well below the estimated accretion rate in Cyg X-3 of 
M acc «5.7x 10- 7 Mq yr" 1 . 
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